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ABSTRACT 



Aims. The properties of the Galactic supernova remnant (SNR) SN 1006 are theoretically reanalysed. 

Methods. Nonlinear kinetic theory is used to determine the acceleration efficiency of cosmic rays (CRs) in the supernova remnant SN 1006. 
The known range of astronomical parameters and the existing measurements of nonthermal emission are examined in order to define the values 
of the relevant physical parameters which determine the CR acceleration efficiency. 

Results. It is shown that the parameter values - proton injection rate, electron to proton ratio and downstream magnetic field strength - 
are determined with the appropriate accuracy. In particular also the observed azimuthal variations in the 7-ray morphology agree with the 
theoretical expectation. These parameter values, together with the reduction of the 7-ray flux relative to a spherically symmetric acceleration 
geometry, allow a good fit to the existing data, including the recently detected TeV emission by H.E.S.S. 

Conclusions. SN 1006 represents the first example where a high efficiency of nuclear CR production, required for the Galactic CR sources, is 
consistently established. 
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1. Introduction 

Cosmic rays (CRs) below the energy 10^^ eV are believed to 
be accelerated i n the shell-type supernov a remnants (SNRs) of 
our Galaxy (e.^. lBerezhko & VoLkl2007l) by means of the diffu- 
sive shock acceleration process. However, this proposition has 
still only a limited observational/theoretical basis. In addition 
SNRs may not only act individually but possibly also more col- 
lectively in the environment of interacting s tellar winds and su 



pernova explosion s in OB associations (e.g. lParizot et al .12004 
Binns et al.|l2007l) . To ensure that Galactic SNRs are indeed 



efficient sources of the Galactic CRs to begin with a number 
of SNRs is required with clearly determined astronomical pa- 
rameters, like the type of supernova explosion, the SNR age, 
the distance, and the properties of the circumstellar medium. 
Applying to such SNRs an appropriate model which consis- 
tently describes the dynamics of the SNR, one can then predict 
the properties of the accelerated particles and the dynamical 
and radiative effects which they produce, like the shock modi- 
fication and the multi-wavelength nonthermal emission. 

The success of the theoretical model can be judged by 
comparison with the experimentally determined overall broad 



band spectrum and with the morphological characteristics of 
the SNR, like its filament structures and the general radial and 
azimuthal variations of the emission pattern, as well as the in- 
ternal dynamics characterized by the contact discontinuity be- 
tween ejected and swept-up mass. From the point of view of 
Galactic CR origin, the key quantity is the efficiency of CR 
production. 

In practice such a program is most often hampered by 
the limited amount and detail of relevant observations. First 
of all, the astronomical parameters of SNRs are as a rule 
poorly known. Even though the SNR age is known for sev- 
eral historical SNRs, the distance is usually quite uncertain. In 
this regard SN 1006 is an exception: the distance was deter- 
mine d using optical mea surements with relatively high preci- 
sion dWinkler et al.ll2003b . 

The second problem is that even the presently most ad- 
vanced nonlinear kine ti c theory of CR accelerat ion in SNRs 



dBerezhko et al.1 Il996t iBerezhko & VoUd 1 19971) Q contains 
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' These papers contain the basic equations used, in the form of a 
transport equation for energetic charged particles in a shock environ- 
ment and its coupling to the gas dynamical equations for the thermal 
plasma. Considering seperately the energetic electron component, a 
corresponding transport equation is added that includes radiative en- 
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physical parameters which can not yet be theoretica lly calcu- 
lated with the necessary precision (see IVoLk:ll2004 for a re- 
view). This concerns the magnitude and the spatial distribution 
of the injection rates of ions and of electrons into the diffusive 
shock acceleration process as well as the extent of magnetic 
field amplification in this process. Fortunately, the values of 
these parameters can be inferred from the observed radio and 
X-ray s ynchrotrori spectra if th ey are measured in sufficient de- 
tail (see lBerezhkoll2005 , 200 81, for reviews). In such a case this 
theory provides a consistent if still approximate model of both 
SNR and CR dynamics, and of the properties of the emission 
produced by the accelerated particles. In particular, the theory 
predicts the high-energy 7-ray spectrum. 

Up to now SN 1006 is the only example, for which 
all astronomical parameters are quite well known (e.g. 
Cassam-Chenai' et a l. 2008). In addition the nonthermal emis- 



injection/acceleration region. Such a strong magnetic field is 
consistent with theoretic al expectation and was confirmed by 



the Chandra detection ( Bamba et al.l 120031 : iLong et alj 120031 



sion in the radio and X-ray ban ds has by now been mea sured 
quite accurately with Chandra jAUen et all |2004 l2008h and 
Suzaku (Bamba_et al. 2008). Beyond that, the VHE emission of 
SN 1006 has been recently detect ed with H.E.S.S., both regard - 
ing its flux and its morphology ("n a umann-Godo et alJ I l2009h . 
This makes SN 1006 uniquely suitable for theoretical study and 
for a detailed comparison with the experimental data. 

We have already applied this nonlinear theory to the case 



of SN 1006 before (Berezhko et al.. 200Z , 


20()3al: Volk et al. 


2003; Ksenofontov et alj 2005; IVolk et al. 


20051 2008bh. In 



the following a brief review of this previous work is given. 
The corresponding model assumes spherical symmetry and a 
quasi-parallel shock. It was, however, argued that efficient nu- 
clear particle injection/acceleration should arise only within 
two polar cap regions, where the SN shock is quasi-parallel 
with respect to the ambient ISM magnetic field w hich in turn 
is approximately perpendicular to the line of sight J Volk et aT 



f2003.> . The calculated size of the efficient nuclear CR produc- 
tion regions, which amounts to about 20% of the shock sur- 
face, corresponds very well to the observed sizes of the bright 
X-ray synchrotron emission regions. Over the remaining 80% 
of the surface the shock is quasi-perpendicular and this leads to 
a depression of the injection of nuclear particles into the diffu- 
sive shock acceleration process. This means that efficient nu- 
clear CR production - that occurs in a spherically symmetric 
model uniformly across the whole shock surface - in the actual 
SNR takes place only within the polar regions. It is therefore 
consistent with a correction for the spherically symmetric so- 
lution for the nuclear CR production and the associated 7-ray 
emission by a renormalization factor /rc ~ 0.2 (for details, 
see the Appendix; to emphasize the three-dimensional charac- 
ter of the configuration a description in terms of magnetic flux 
tubes is given.). In the first study (Bere zhko et al. 2002) it was 
demonstrated that the observed properties of the nonthermal 
emission can be consistently understood assuming a consider- 
ably amplified magnetic field Bd ~ 120 fiG within the nuclear 



ergy losses jBerezhko et al.ll2002h . Also the particles' diffusion prop- 
erties in the form of the Bohm limit, wherever nuclear particles are 
effectively injected into strong shocks, as well as the assumed heating 
rate of the thermal plasma by wave dissipation are detailed in these 
papers. 



of filament structures in th e nonthermal X-ray em ission of the 
rims of the SNR shell (see iBerezhko et al.ll2003ai for an inter- 
pretation). A filamentary structure of this character is indeed 
expected in young SNRs, where the magnetic field is strongly 
amplified, as the result of the strong synchrotron losses of the 
X-ray emitting CR electrons. In fact, since there exists gener- 
ally a simple relation between the filament thickness and the 
interior magnetic field strength dBerezhko & Volkll2004al) . the 
measurement of the filament thickness represents a new and in- 
dependent method for the determination of the magnetic field 
strength inside young SNRs. It is important to note that the 
interior magnetic field strength, determined in such a way for 
all known Galactic SNRs, agrees very well with its value de- 
rived from a fit of the shape of the overall synchrotron spec - 
tra, wherever both methods can b e appl ied dVolk et alj|2005l) . 
It was also shown jBerezhko et al ] l2002h that nonlinear kinetic 
theory of CR acceleration is consistent with the TeV-emission 
detected by CANGAROO in 1998 for a value A^h = 0.3 cm'^ 
of the ambient interstellar medium ISM density from the range 
0.05 < A^H < 0.3 cm~^ existing in the literature. However, 
SN 1006 was not confirmed as a TeV source by the H.E.S.S. 
experiment in a total of 18. 2h (in 2003) and 6.3h (in 2004) 
"livetime" of ON source observat ions with SN 1006 in the 
field of view ( Aharonian et alj|2005i) . The H.E.S.S. upper limit 
was roughly one order of magnitude lower than the published 
CAN GAROO flux. In a subseq uent paper by the present au- 
thors (iKsenofontov et al" I l2005 i) it was demonstrated that this 
H.E.S.S. upper limit does not invalidate the theoretical picture 
on which the previous calculation of the 7-ray emission spec- 
trum had been based. As it was shown, it is rather the value 
of the external astronomical parameter iVn that strongly influ- 
ences the hadronic 7-ray flux: the hadronic 7-ray flux is very 
sensitive to the ambient gas density iVn and the H.E.S.S. up- 
per limit requires A^h < 0.1 cm^'^. It is important to realize 
that only a rather strongly amplified magnetic field, produced 
nonlinearly by efficiently accelerated nuclear CRs, is compat- 
ible with such a low 7-ray flux which is therefore expected to 
be predominantly of hadronic origin. Otherwise the overpro- 
duction of accelerated electrons would leed to an unacceptably 



high inverse Compton (IC) 7-ray emission (IBerezhko et al 



2002h . This requires a strong component of accelerated nuclear 
particles whose energy density substantially exceeds that of 
the synchrotron electrons. The opposite simplest IC scenario, 
which is based on the assumption that the entire nonthermal 
SNR emission is produced by accelerated electrons without any 
substantial nuclear CR production and magnetic field amplifi- 
cation, substantially overpredicts the exi sting 7-ray up per lim- 
its for aU known type la Galactic SNRs dVolk et al.lEo08b.) . 

Sin ce TeV 7-ray emission of SN 1006 was recently de- 
tected dNaumann-Godo et al. I I2OO9I) it is worthwhile to recon- 
sider SN 1006 on the basis of thi s nonlinear theory. Comp ared 
with the previous consideration (IKsenofontov et al.ll2005h this 
makes it possible to determine the ambient ISM density within 
narrow limits. In addition, we shall demonstrate that the values 
of the relevant physical parameters - proton injection rate, elec- 
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tron to proton ratio and downstream magnetic field strength - 
are determined for SN 1006 with appropriate accuracy from the 
measured synchrotron spectrum. The 7-ray morphology agrees 
with the theoretical expectations regarding the morphology of 
ion injection and t he corresponding morphology of magnetic 
field amplification (IVolketalJ2003l . It is therefore indeed con- 
sistent with a correction for the total nuclear particle pressure 
and the total 7-ray flux obtained in the spherically symmetric 
solution by a renormalization factor ~ 0-2. The coiTec- 
tion implies the full value of particle pressure and 7-ray flux, 
as calculated in the spherically symmetric model, in the quasi- 
parallel part of the magnetic flux tubes that thread the polar cap 
regions, but a value of zero elsewhere. The above parameter 
values plus the mentioned morphological aspects allow a physi- 
cal explanation of all existing multi-wavelength data, including 
the H.E.S.S. measurement of the TeV spectrum. There is pos- 
sibly one exception to this conclusion. It concerns the radius of 
the contact discontinuity between ejected and swept-up mass 
relative to the ra dius of the SNR blast w ave, cf. recent data dis- 
cussed by ICassam-Chenai' et al.l (l2008h which indicate a large 
value - apparently quite clos e to unity - o f this ratio. However 
another very recent analysis (iMiceli et al ]'2009) claims values 
^ 0.9 for this ratio which are rather in line with quantitative 
expectations for a CR-modified shock in a n object of compara 



ble evolutionary phase like Tycho's SNR dWarren et al.ll2005 



Cassam-Chenai et al. 2007: Volk et al. 



2008al) . 



Future work 



will have to assess these differences. 

2. Results and discussion 

Since SN 1006 is a type la supernova it presumably expands 
into a uniform ISM, ejecting roughly a Chandrasekhar mass 
Moj = 1.4M0. Since the gas dens ity indeed varies only mildly 
across the SNR ( Acero et al.l2007 ). it appears reasonable to as- 
sume also the circumstellar magnetic field to be uniform. The 
ISM mass density po = l-4mpAfH, which is usually character- 
ized by the hydrogen number density TVh, is an important pa- 
rameter which strongly influences the expected SNR dynamics 
and nonthermal emission. 



As in our earlier study (IKsenofontov et alJl2005l) we solve 
here the coupled set of nonlinear equations, mentioned before, 
in order to find the optimum set of physical parameters of 
SN 1006 which gives a consistent description of the observed 
overall dynamics and of the nonthermal emission together with 
its morphology. The theory includes all the important physical 
factors which influence CR acceleration and SNR dynamics: 
shock modification by CR backreaction, MHD wave damping 
within the shock transition, a consistently determined CR spec- 
trum, and the spatial distributions in each evolutionary phase. 
In addition it includes synchrotron losses of CR electrons and a 
determination of all nonthermal emission processes, produced 
in SNRs by the accelerated CRs. It had also been shown that the 
values of the key parameters (proton injection rate 77, electron 
to proton ratio Kcp and interior (downstream) magnetic field 
strength B^) which can not be predicted theoretically with the 
required accuracy, can be determined from a fit of the observed 
synchrotron emission data. It is of basic importance here that 
the parameter values for SN 1006, determined in this way, were 



very well confirmed by the Chandra meas urements of the fine 
structure of the n onthermal X-ray emission dBamba et alJ2003l: 

Long et al ] |2003) . as analyzed by lBerezhko e t al. ('2003aV 

Compared with the previous study fesenofontov et alj 



2005h th e analysis includes the most accur ate X-ray data from 
Chandra (" Allen et al.l2004 ') and Suzaku gv feamba et al.l2008l) 
that make it possible to quite precisely determine the accept- 
able range of these parameter values. Also a time-dependent, 
amplified upstream magnetic field 



Boit) = Boitsn)VPcit)/Pc{tsn) 



(1) 



is used - in the spherically symmetric calculation over the en- 
tire shock surface - where Pc is the CR pressure and tsn is the 
SNR age (see .Berezhko & V61k 2004b) . The value of the mag- 
netic field at the current epoch Bq {tsn) is determined as a result 
of the best fit of the measured synchrotron spectrum (see be- 
low). In the downstream region the magnetic field is assumed to 
be frozen into the gas with the postshock field value B2 ~ ctBq, 
where a is the shock compression ratio. Under these assump- 
tions the magnetic field strength varies only slightly within 
the downstream region of swept-up ISM. Therefore below the 
value _Bd = B2 is used as the interior magnetic field strength. 
Also the notations B^ and Bq are used for their values at the 
current SNR epoch t ^ tgn- 

Since the properties of accelerated CR nuclear and electron 
spectra and their dependence on the relevent physical parame- 
ters, as well as the dynamical properties of the system were de- 



scribed in detail in previous papers (e.g lBerezhko et al.ll2002 



Ksenofontov et al.ll2005) . they will not be discussed here (see 
however the Appendix for a summary). 

The values of the SN explosion energy E^n = 1.8 x 10^^ erg 
and Esn — 1-5 x lO*^^ erg a re taken to fit the o bserved shock 
size Rs and shock speed Vs (IMoffett et al.lll993h at the current 
epoch t 10^ yr (see Fig.l) for the ISM hydrogen number 
densities = 0.05 cm'^ and iVn = 0.035 cm respec- 
tively. These densities are consistent with the observed level of 
the VHE emission (see below). Note that the calculation pre- 
sented in Fig. 1 corresponds to the best fit value of the amplified 
upstream magnetic field strength Bq = 30 /iG, even though 
the shock size Rg and shock speed Vs are quite insensitive to 
Bo. The resulting cuiTent total shock compression ratios a for 
iVn = 0.05 cm-3 and = 0.035 cm^^are a = 4.9 and 4.7, 
respectively, whereas the subshock compression ratios CTs are 
both close to as = 3.7. 

As the most rel iable estimate for the distance the value d = 
2.2 kpc was taken (' WinkleretaLll2003h . 

All these results stem of course mainly from the gas dy- 
namics part of the governing equations. 

Fig. 2 shows the volume-integrated momentumg distribu- 
tions of energetic protons (representative for the nuclear com- 
ponent) and electrons. The proton cutoff energy e^^^ is above 
10^^ eV. The electron spectral index 7 « 2 equals that of the 
protons below the onset of synchrotron losses at the energy 
e; « 1 TeV. The electron cutoff energy e^^^ i^^^" 10^^ 
Figs. |3] Eland |5] illustrate the consistency of the synchrotron 
spectrum, calculated for the ISM number density A^h = 
0.05 cm~^ with the best set of parameters (77 = 3 x lO^"', 
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Fig. 1. Gas dynamical evolution for the two values of the total 
explosion energy Egn = 1-8 x 10^^ erg and Esn — 1-5 x 10^^ 
that correspond to the values A^h = 0.05 cm~'^ and A^h = 
0.035 cm"'^ of the ISM hydrogen number density, respectively. 
The quantities Rs and Vs denote the shock radius and the shock 
velocity, their time evolution being fitted to the observed values 
shown. The vertical dotted line indicates the present epoch. 




Fig. 3. Differential synchrotron radio emission flux as a func- 
tion of frequency, calculated for the ISM hydrogen number 
density A^h = 0.05 cm""^ at the proton injection rate 77 = 
3 X 10"'', electron to proton ratio, Kcp = 4 x lO^"*, and mag- 
netic field strength Bq = 30 upstream of the sub shock . The 
observed radio emission flux from a com pilation of lAIlen et al 
(|2004 (see also lAUen et al]|200lll2008l) is shown. 




log(p/mpC) 

Fig. 2. Overall (=spatially integrated) spectra of energetic parti- 
cles. Solid and dashed lines correspond to protons (p) and elec- 
trons (e), respectively. 



Kcp = 4 X 10^'', Bo = 30 fiG), with the observed spatiaUy in- 
tegrated spectra. Note that we use the differential synchrotron 
fluxes in the radio (Su) and in the X-ray (dF^/de-y) range in a 
slightly different form. They are however simply connected by 
the relation vS^ ^ e^dF-y/de^. 

It should be noted that the values a > 0.5 of the differen- 
tial radio spectral index a — —dlnSi./dlni', as observed in 
young SNRs, require a curved electron spectrum that hardens 



to higher energies as predicted by nonline ar shock accelera- 
tion m odels, as originally emphasized by Revnolds & Ellison 
(Il992r) . To have a — 0.57 in the radio range, as observed 
for SN 1006, one needs efficient CR acceleration with a pro- 
ton injection rate 77 = 3 x 10^''^ which leads to the required 
shock modification, a nd a high interior magnetic field i?d > 



130 mG is req uired (Be rezhko et al. 2002; 



Volk et al. 2005 



Ksenofontov et a l. 2005). The hardening of the observed radio 



synchrotr on spectrum for S N 1006 has been recently demon 
strated bv^ AUen et al. ( 2008 ). Detailed X-ray synchrotron spec 
tral measurements are, however, required t o find the optimum 



value of the magnetic field strength B^ dKsenofontov et al 



2005h since for a given fit of the calculated synchrotron spec- 



trum in the radio range the X-ray synchrotron amplitude is very 
sensitive to Bd. The reason is the foUowing: in the polar re- 
gions, where the scattering mean free path decreases down to 
the particle gyroradius (Bohm limit), the spectrum of acceler- 
ated electrons at the shock front is quite insensitive to the mag- 
netic field strength B^ if all other relevant parameters are fixed; 

— 1/2 

the exception is their maximal ener gy e^„^ oc VsB^ which 
is due to synchrotron losses (e.g. iBerezhko & Volk l2004bh . 



Since all energy of the electrons at high energies e > ei oc B^^, 
where synchrotron losses are significant, is rapidly and com- 
pletely transformed into the synchrotron emission, the same is 
true for the high frequency (i' > iy\ oc ef / B^) part of the syn- 
chrotron spectrum including its cutoff part, because the cutoff 
1/2 

frequency Vms^x oc efnax^d doss not depend on B^- (In the 
case Bo ~ 30 /iG and B^ = 150 /iG we have e^jax ~ 10^^ eV 
and i^max ~ 10^* Hz.) However, varying the magnetic field 
strength B^ changes the amplitude of the overall calculated 
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Fig. 4. X-ray synchrotron emission flux as a function of photon energy for the same case as in Fig. 1 (solid line). Calculations for 
i?o — 27 /iG (dashed line) and for Bq — 33 jjG {dash-dotted line) are also presented. Th e observed X-ray flux from Chandra in 
red color for a small region of the bright northeas tern rim of SN 1006 (lAllen et al.ll2004l) . and the global X-ray spectrum in blue 
color as observed by Suzaku feamb a et al. 2008 ) are shown. The Chandra X-ray flux was normalized such as to be consistent 
with the Suzaku flux for energies e > 2 keV. 



synchrotron spectrum at high frequencies > ^-i in order to 
keep the fit of the measured radio spectrum at low frequencies. 
Therefore the high frequency part of the spectrum changes 
due to the variation of the breaking frequency v\ cx B^^ . Since 
below v\ the synchrotron spectrum has the form Si, cx i^^" with 
a ~ 0.5 we conclude that the high frequency part of the syn- 
chrotron spectrum depends on Bd like Su{v\) oc B^ The 
amplitude Sy{v\) is therefore the main determinant of the value 
of i?d- Together with the observed synchrotron amplitude in 
the radio frequency region, where radiative losses play no role 
yet, also the amplitude of the nonthermal electron distribution 
and thus the electron/proton ratio ratio Kf,-p is then determined. 
Note tiiat the synchro tron spectrum is almost i nsensitive to the 
ISM density TVr (e.g.' Ksenofontov et alj|2005l) . 

To obtain the fit to the data shown in Fig. [3] and Fig. |4] 
it is asserted that only the highest-energy part of the ob- 
served global X-ray spectrum, corresponding to e > 2 keV, 
is of predominantly nonthermal origin. Towards lower ener- 
gies e < 2 keV the contribution of thermal emission to this 



X-ray spectrum progressively increases, as clearly indicated 
by the Suzaku spectrum. For these considerations the X-ray 
Chandra flux fr om a small region of the bright north eastern 
rim of SN 1006 dAUen et al.ll2004 is used, where t he contribu- 



tions from ther mal X-rays are presumably minimal (lAllen et al 
2004; Cassa m-Chenai' et al. 2008). In order to apply it to the 
whole remnant, this X-ray flux was normalized to be consis- 
tent at energies e > 2 keV w ith the global X-ray spectrum as 
observed by Suzaku (IBamba e t al. 2008). This consistency in 
fact exists. For quantitative purposes the normalised Chandra 
spectrum above 2 keV is used, under the assumption that it 
corresponds to pure nonthermal emission. Comparison with the 
experimental X-ray data shows in Figs. [3] and H] that the opti- 
mum magnetic field value is So = 30 /iG, and corresponds to a 
downstream field B^ ~ 150 /iG. This is in agreement with the 
field amplificat ion that is implied by the filamentary structures 
in hard X-rays dVolk et al.ll2005h . 

Fig- 13 shows the contours of characterizing the quality 
of fit of the radio spectrum, and of the X-ray spectrum for e > 
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Fig. 5. Contours of equal deviations Ax^ (shown on the curves) 
from the minimum value of ^ , in the plane that is spanned by 
the electron to proton ratio ifop and the proton injection rate r\. 
The asterisk corresponds to the best fit values 77 = 3 x 10^^ 
and JsTcp ~ 4 x 10^* with the minimum value x^/dof =1.3. 
The fit uses jointly the radio and nonthermal X-ray data. 



2 keV, for different values of 77 and K^^. One can see that the 
quality of the fit is rather good, allowing only a rather small 
range of these parameters. 

The quality of fit of the radio spectrum, and of the X-ray 
spectra for e > 2 keV, for Bq = 30 is characterised by the 
value x^/dof = 1.3. An increase Ax^ « 1 of x^ (see Fig. HI 
implies a change of the magnetic field strength value by only 
10%. 

The 7-ray m orphology, as found in the H.E.S.S. mea- 
surements (.Naumann-Godo et al. I l2009l) . is consistent with the 
prediction of a polar cap geometry by IVolk et al.l j2003b . 



Such a geometry has been also found experimentally from 
an analysis of th e synchrot ron morphology in hard X-ra ys by 
Rothenflug eTal] (ll004|) and lCassam-Chenai' et alj (l2008h . This 
means that the 7-ray emission calculated in the spherically 
symmetric mo del must be renormalized (reduced) by a factor 
/lo ~ 0.2, as in lKsenofontov et alj (120051) . The renormalization 
factor is applied here. 

This morphology is also the key argument for the exis- 
tence of an energetically dominant nuclear CR component in 
SN 1006, because only such a component can amplify the mag- 
netic field to the observed degree. If on the contrary the existing 
energetic electron component would have to drive the magnetic 
field amplification all by itself then we would, as a minimum, 
require ^ i?g/(87r), where P^, is the pressure of CR elec- 
trons at the shock front. Clearly P^ ^ K^pPp, where Pp « Pc is 
the pressure of CR protons and Pc is the pressure of the nuclear 
component. Since Pc « O-^poV^ we have Pc « 0.07P^/(87r) 
for Na = 0.05 cm^^ and T4 w 4000 km/s. Therefore ac- 



celerated electrons are not able to amplify the magnetic field 
to the required level by a large margin. On the other hand, 
the effectively accelerated nuclear component has the pressure 
Pc ~ 200Po/(87r) and can therefore readily amplify the field 
from a purely energetic point of view. 

The amplified field must also be the reason for the corre- 
sponding polar cap-type morphology of the synchrotron emis- 
sion in hard X-rays. Even though the weak radio synchrotron 
e mission, essen tially all around the periphery of SN 1006 



dRothenflug et a l. 2004: iCassam-Chenai' et al.l l2008h . demon 



strates that electrons are at least to some extent injected ev- 
erywhere into the acceleration process over the shock, they 
reach the multi-TeV energies for X-ray synchrotron emission 
only in the polar caps, where the amplified field allows their 
acceleration to these energies. In the equatorial region the up- 
stream magnetic field strength corresponds to its ISM value, 
at the height of SN 1006 above the Galactic plane probably 
not more than 3/iG. Particle scattering is comparatively weak 
and the maximum electron energy is expected to lie below the 
limit, where the acceleration rate equals the synchrotron loss 
rate. Then the maximum electron energy is - like that of the nu- 
clear particles - determined by the system g eometry jBe rezhki 
1996); an explicit calculation can be found in lBerezhko & Vol 
(i2004 b). Assuming then for simplicity that the spatial diffusion 
coefficient has the form k — S ■ KBohm, with a constant factor 
S > 1, requires ^ > 4 even for this low field strength in or- 
der to reduce the equatorial maximum synchrotron frequency 
t'max to t'niax < 10^^ Hz, for an azimuthally uniform shock 
velocity 14 ~ 4000 km/s. Such a plausible reduct i on of i^max 
is c onsistent with the fi ndings of iRothenflug et al.l (120041) . see 
also lMiceh eTall (|2009|) . 

The only important parameter which can not be determined 
from the analysis of the synchrotron emission data is the exter- 
nal density iVn. Therefore we have performed the calculations 
for the pair of values iVn = 0.05 cm^'^ and iVn — 0.035 cm~'^ 
which seems to bracket the density range consistent with the 
H.E.S.S. 7-ray measurements. 

Fig. |6] shows the total (7r°-decay plus inverse Compton 
(IC)), and seperately the 7r°-decay and the IC 7-ray energy 
spectra of the remnant, calculated for iVn = 0.05 cm^^^ and 
Nu = 0.035 cm-3. 

H.E.S.S. has reported the measured value of the energy flux 
^{^1,^2) = 2.5 X 10^^^ erg/(cm^ s) of 7-rays with energies 
ei < £7 < £2, where ei = 0.2 TeV and £3 = 40 TeV. It is 
related with the differential 7-ray energy spectrum according 
to the expression 



«'(ei,e2) 



dF 



In 



£2 



< e 



2 '^-^7 



> 



In a selfconsistent picture the theoretically predicted form 
of the differential spectrum dFj/de^, cf. Fig. 6 is used to cal- 
culate the mean experimental value of the differential energy 
flux < e^^dFj/dcj >— <i>(ei, e2)/ln(e2/ei), shown in Fig.|6] 

Note that the peak in the 7-ray energy spectrum at e-^ « 
5 TeV and the quasi-exponential cutoff at higher energies is 
due to the substantial contribution of the IC component: at 
TeV energies the IC component contributes about half of the 
total 7-ray flux for Nn — 0.05 cm^^, whereas for A'h = 
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Fig. 6. Total (7r°-decay + IC) {black lines), 7r°-decay (blue 
lines), and IC (green lines) differential 7-ray energy fluxes as a 
function of 7-ray energy, calculated for ISM hydrogen number 
densities A'h — 0.05 cm~'^ (solid lines) and Nn = 0.035 cm^-^ 
(dashed lines) for the parameters rj — 2.9 x 10^"*, Kep — 
4.1 X 10~^ and Bd = 150 ^G, derived from the fit of the syn- 
chro tron spectrum. The H.E.S.S. value dNaumann-Godo et al 
l2009i) is shown as well. 



0.035 cm^'^ it is already two thirds. Since the maximal pro- 
ton energy, reached at some time during the evolution, is about 
Cmax ~ 2 X 10^^ eV, the corresponding 7r°-decay component 
has a power law form dF^/de^ oc e~'^ with 7 2 up to the 
cutoff energy e-, « 10^^ eV. The cutoff energy is defined here 
as that energy, where the spectral energy density has dropped 
by a factor 1/e from its maximum value which given by log 



eldF-y/de-^ 



-0.8. 



According to Fig. |6] the H.E.S.S. data are consistent with 
an ISM number density from quite a narrow interval 0.035 < 
A'h < 0.05 cm^'^, since for the theoretically derived 7- 
ray spectrum we have $ = 2.1 x 10^^^ erg/(cm'^ s) and 
$ = 2.9 X 10^12 erg/(cm2 s) for Nu = 0.035 cm^^ and 
TVh = 0.05 cm~'^ respectively. It should be noted that the cor- 
responding explosion energy Egn ~ 1-7 x 10^^ erg is close to 
the upper end of the typical range of type la SN explosion en 



ergies that vary by a factor of about two dGamezo et al.ll2005 



Bhnnikov et al 



Aceroetal 



20061) 



(120071) find the value Nu « 0.05 cm"^ on the 
basis of X-ray measurements. The above interval is consistent 
with their result. 

From Fig. 6 the 7-ray spectrum produced by the nuclear 
CRs is very close to the IC emission spectrum produced by 
CR electrons alone. Since the differential energy spectrum of 
freshly accelerated nuclear particles and electrons is rather 
close to a spectrum TVe oc e~^, and since the electrons with 
energies e > ei ^ 1 TeV subsequently undergo significant 
synchrotron cooling in the interior, leading to the spectrum 



c-3 



not only the amplitude but also the shape of these 



two components are very similar within the energy interval 
10^° < e-y < 10^^ eV. Therefore the VHE 7-ray spectrum 



alone is not able to discriminate between the hadronic 7r°- 
decay and the lept onic IC 7-rav components . However, it was 
already shown by ' Ksenofontov et al] ( 2005 ) that such a low 
VHE emission flux, with a highly depressed IC 7-ray flux, is 
only possible if the nuclear CR component is efficiently pro- 
duced with accompanying strong magnetic field amplification. 
In the framework of the interpretation developed in this pa- 
per the most direct evidence for the energetic dominance of 
a nuclear energetic particle component is the observed 7-ray 
morphology. It corresponds to the theoretical prediction and is 
consistent with all other measurements. 

A last point concerns the radial extent i?c of the contact dis- 
continuity between ejected and swept-up mass relative to the 
radius Rs of the SNR blast wave, cf. rece nt data presented and 
discussed by Cassam-Chenai' et al. ( 200 8^. The ratio Rs/Rc is 
given as 1.04 ± 0.03 outside the synchrotron rims, and essen- 
tially as 1 in the region within the synchrotron rims. In the case 
of Tyc ho's SNR we have discussed in quantitative detail (e.g. 
IVolk et al. 2008a) the reduction of this ratio compared to a pure 
gas shock as a result of the considerable shock modification 
produced by accelerated nuclear CRs, which leads to the in- 
crease of the shock compression. Qualitatively such considera- 
tions agree with the experiment for SN 1006: the ratio Rs/Rc is 
larger in the equatorial region, where CR injection/acceleration 
is inefficient, and it is essentially smaller within the polar re- 
gion, where CRs are efficiently produced. While for Tycho's 
SNR particle acceleration gives a good theoretical explanation 
of the relatively small ratio within our model, this is clearly 
impossibl e quantitatively for the above numbers in the case of 
SN 1006. Cassam-Chenai' et all (l2008b believe that the value 
of the contact discontinuity radius R^ was somehow overesti- 
mated "... since our measurements are likely to be affected by 
projec tion and other effec ts, On the other hand, very re- 
cently ^iceiretaL ( 20091) have found numbers Rs/Rc of the 
order of 1 . 1 which are in the expected range. Future work will 
have to resolve this difference. 

3. Summary 

Since the relevant astronomical parameters as well as the syn- 
chrotron spectrum of SN 1006 are measured in impressive de- 
tail it is possible to determine the values of the relevant physical 
parameters with the appropriate accuracy for this SNR: proton 
injection rate rj = (2.9 ± 0.6) x 10~^, electron to proton ra- 
tio Kep = (4.1 ± 0.3) X 10~^ and downstream magnetic field 
strength Bd = (150 ± 15) ^iG. 

As a result the flux of TeV emission detected by H.E.S.S. 
is consistent with the ISM number density 0.035 < iVn < 
0.05 cm~^. The corresponding hydrodynamic SN explosion 
energy Esn = 1.7 x 10^^ erg is close to the upper end Esn = 
1.6 X 10^^ erg of the typical range of type la SN explosion 
energies that vary by a factor of about two. Also the magnetic 
field amplification properties of this SNR are well understand- 
able as the result of azimuthal variations of ion injection over 
the projected SNR circumference and corresponding accelera- 
tion which lead to a polar cap-type morphology for the X-ray 
synchrotron as well as the 7-ray emission. As a consequence, 
the recent H.E.S.S. measurement of a dipolar morphology also 
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of the TeV 7-ray emission is considered as the most significant 
empirical evidence for the existence of an energetically domi- 
nant nonthermal, nuclear component. 

In conclusion SN 1006 appears to represent the first exam- 
ple where a high efficiency of nuclear CR production, required 
for the Galactic CR sources, is consistently established. 
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Appendix A: 

The relation between the model calculations in spherical sym- 
metry (ID) and the physical characteristics of the system in 
the presence of a symmetry-breaking magnetic field in the 



circumstellar medium have been discussed earlier (iVolk et al 



2003h . In order to clarify the necessary corrections to the ID- 
calculation for SN 1006 a qualitative discussion of the physical 
situation is given here. 

A key feature is the empirical fact that the force den- 
sity of the magnetic field in the diffusive shock acceleration 
process can be neglected compared to the gas ram pressure, 
the thermal pressure and the CR pressure. The average mag- 
netic field configuration can thus in principle be calculated 
purely kinematically from the induction equation. However, 
the magnetic field direction influences the injection at least 
of heavy ions in a significant way. In a quasi-parallel shock, 
where the angle QnB between the locally mean field and the 
shock normal is small compared to 90°, the required velocity 
of suprathermal ions from the downstream region, required to 
be able to cross the shock into the upstream region, is much 
smaller than in a quasi-perpendicular shock, where Q^b is 
close to 90°. For a steeply falling velocity distribution of the 
shocked downstream gas this leads to effec tive suppression of 
ion injection in quasi-p erpendicular shocks ( Ellison et al.ll995 ; 



Malkov& VoIkllI995l) . The acceleration of injected particles, 
on the other hand is efficient for all 8nB for non-relativistic 
shocks, as long as cos9nB ^ K/c, where Vs and c denote 
the sh ock velocity and the speed of light, respectively (IDrury 
19831). The injection of electrons is not yet well understood. 
From the observed azimuthal di stribution of the radio syn- 
chrotron emission in SN 1006 (e.g. lCassam-Chenai et aljioosi) 
electron injection seems to operate for all angles Qnh, although 
apparently less efficiently in the quasi-perpendicular region. 
This is the least restrictive assumption and it is made here in 
the sequel. 

A second point is that in the instantaneously quasi-parallel 
part of the magnetic flux tubes, in which ions are injected, 
the accelerated particles drive a current against the flow of 
the upstream gas that generates strong magnetic fluctuations 
and amplifies the magnetic field, as refeiTed to in the main 
text. This enhances the acceleration and leads to a pres- 
sure Pc of the energetic nuclear particles that is compara- 



ble to the gas ram pressure for strong shocks, e.g. in young 
SNRs, as well as to heating of the gas by wave dissipation, 
both of the non-resonant and of the resonant Alfvenic modes 
(McKenzie & Volk 1982; Lucek & Bell 2000; Bell & Lucelg 
[200ll:lBeli2004l:IPeUetier et aI.ll2006HzirakashviIi et all2008l) . 
Especially the formation of density fluctuations and corre- 
sponding secondary sho cks is viewed here as a reason fo r 
strong wave dissipation ( Berezhkol 2008 : Volk et al. 2008ah . 
The pressure of accelerated electrons is negligible in compar- 
ison. This implies a strong hadronic 7-ray emission from that 
part of these flux tubes, which became at some time quasi- 
parallel. This hadronic flux is negligible elsewhere because 
there are few energetic nuclei there. A certain amount of cross- 
field diffusion will, however, also populate neighboring flux 
tubes and somewhat smear out the acceleration boundaries 



dVoIk et al.ll2003h . Accelerated electrons will emit strong syn- 
chrotron radiation in the enhanced magnetic field regions, but 
may also yield a more moderate synchrotron flux elsewhere. 
The reason is that by assumption they are injected and ac- 
celerated to some extent "everywhere" on the shock surface, 
albeit to considerably lower maximum energies in the quasi- 
perpendicular regions. This is the result of the lack of self- 
excited magnetic field fluctuations and of field amplification 
(see section 2). 

The third point regards differences in overall dynam- 
ics of the quasi-parallel and quasi-perpendicular regions. 
Conservation of the overall mass, momentum and energy fluxes 
allows only a different partition within the momentum and en- 
ergy fluxes between their kinetic, thermal and nonthermal com- 
ponents, but the sums of these fluxes always remain the same. 
However, the compressibility of the relativistic CR "gas" is 
higher than that of the nonrelativistic thermal gas, and therefore 
this different partition implies dynamical differences between 
the quasi-parallel and the quasi-perpendicular regions. 

When primarily the cumulative effect of the higher com- 
pressibility counts, like in the relative distance between the 
outer SNR blast wave and the following contact discontinuity, 
then the different partitions in a CR-modified shock and a pure 
gas shock play the determining role. 

However, in the present context, where the main character- 
istics of the radiation effects and of their spatial distributions 
are to be evaluated, such dynamical differences can to first ap- 
proximation be neglected. In practice, for SN 1006, the over- 
all shock compression ratio changes from a = 4 without CR 
production to ct ~ 5 with efficient CR production, given the 
assumed strong gas heating in the shock precursor 

Thus, for the evaluation of the hadronic 7-ray emission, 
only those flux tubes from the ID-configuration gas dynamic 
configuration are considered which connect to the part of the 
present phase shock surface, where ion injection is efficient. 
This requires the correction factor f^c < 1. Since the particle 
transport along magnetic field lines is primarily diffusive, adi- 
abatic expansion of the thermal gas leads to adiabatic losses 
of the energetic particles and weakens the interior particle dis- 
tribution. In addition almost all flux tubes have been orig- 
inally quasi-perpendicular. The coiTesponding parts of these 
flux tubes are now in the deeper interior of the SNR and there- 
fore, to first approximation, lack there accelerated nuclei. For 
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electrons also the radiative losses increase towards the interior 
where the 'older' particles reside. These arguments show that 
both the synchrotron and the 7-ray emission are concentrated 
near the shock surface rather than in the interior As long as 
the ejected mass is not negligible compared to the swept-up 
mass, the ejecta volume in the interior with its extremely small 
magnetic field strength will appear as a region with particularly 
little radiation. For SN 1006 with its supposedly simple mag- 
netic field geometry, the evaluation of the ion injection fraction 
/re of the shock surface has been given by Volk et aL (2003). 

Also the electro n spectrum is ca l culated in spherical sym- 
metry from Eq.3 of B erezhko et al. I (12002 '). However the ac- 
tual evolutionary phase of SN 1006 appears to be still rather in 
the transition from the sweep-up phase to what one can losely 
call a Sedov phase. Being pushed away by the ejecta the den- 
sity of the shocked circumstellar medium has remained rather 
uniform. In the quasi-parallel flux tubes, downstream of the 
shock, with massively injected nuclear particles, the field 
should therefore be strong and roughly uniform if field dissi- 
pation is negligible. In the quasi-perpendicular regions on the 
other hand, the field is only about four times larger than in the 
ambient circumstellar medium. Finally, in the ejecta mass the 
magnetic field strength is presumably very small at the present 
time. Taken all this together, the overall magnetic field strength 
in the SNR is far from uniform. Yet what counts primarily for 
the integrated synchrotron emission are the quasi-parallel re- 
gions. 



It had indeed been assume d in past work dBerezhko et al 



2002tlKsenofontov et al ]l2005h that to first approximation only 



those magnetic flux tubes play a role in the spatially integrated 
emission which end at the forward shock in the polar caps. 
Therefore this integral already contains the correction factor 
/rc implicitely and can be directly compared to the observed 
integral flux. The key consideration in all of this is that the SNR 
consists of a strongly CR-modified part with strong magnetic 
field amplification and a part which is essentially unmodified. 
The correction factor changes only the amplitude but not the 
form of the calculated spectrum. 

Adding a small population of electrons accelerated in the 
unmodified shock of the quasi-perpendicular regions would 
smear out the specifically nonlinear properties of the syn- 
chrotron spectrum. Therefore, the observed synchrotron spec- 
trum exhibits the minimum nonlinear modifications which in 
reality should be even somewhat stronger This means that 
the derived ion injection rate, magnetic field amplification and 
electron/proton ratio are lower hmits for the quasi-parallel re- 
gions. 
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